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^ \ Abstract 

\ We solve the equations of motion for a scalar field with domain wall boundary 

\ conditions in a Friedmann-Robertson- Walker (FRW) spacetime. We find (in agreement 

with Basu and Vilenkin) that no domain wall solutions exist in de Sitter spacetime 

, for h = H/m > 1/2, where H is the Hubble parameter and m is the scalar mass. 

P^, In the general FRW case we develop a systematic perturbative expansion in h to 

Oh! arrive at an approximate solution to the field equations. We calculate the energy 

(D . 

, momentum tensor of the domain wall configuration, and show that the energy density 

can become negative at the core of the defect for some values of the non-minimal 
^ . coupling parameter ^. We develop a translationally invariant theory for fluctuations 

^ . of the wall, obtain the effective Lagrangian for these fluctuations, and quantize them 

using the Bunch-Davies vacuum in the de Sitter case. Unlike previous analyses, we find 
that the fluctuations act as zero-mass (as opposed to tachyonic) modes. This allows 
us to calculate the distortion and the normal-normal correlators for the surface. The 
normal-normal correlator decreases logarithmically with the distance between points 
for large times and distances, indicating that the interface becomes rougher than in 
Minkowski spacetime. 
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1 Introduction 



The formation of topological defects in a cosmological phase transition is a basic result in 
particle cosmology lore H. The existence of defects provides us with a powerful tool 
which both constrains particle theory models, such as those containing domain walls |Q 
and magnetic monopoles 0], and solves some outstanding cosmological problems, such as 
the formation of structure with the use of cosmic strings H. Furthermore, during out- 
of-equilibrium phase transitions, topological defects such as interfaces (domain walls) are 
important ingredients for the dynamics of phase separation and phase coexistence. 

It has always been assumed that if a field theory defined in Minkowski spacetime admits 
topological defects, the same will be true in an expanding universe, and that furthermore, 
there will not be any significant differences in the physical characteristics of these defects. 
Recently, however, Basu and Vilenkin p have shown that defects in a de Sitter background 
can have properties that are quite different from defects in fiat spacetime. In this work, we 
continue and extend this program by analyzing the equations of motion for a domain wall in 
a background Friedmann- Robert son- Walker (FRW) spacetime. We find that although these 
equations cannot be solved analytically in general, a systematic perturbative expansion can 
be set up, with the small parameter h = H/m being the ratio of the correlation length 
of the field to the horizon size. Such a perturbative expansion arises naturally from the 
requirement that the width of the domain wall be much smaller than the particle horizon, 
to allow for kink boundary conditions within the microphysical horizon. The zeroth order 
solution is taken to be the standard flat spacetime domain wall conflguration, except that it 
is a function of the physical spatial coordinate. The higher order terms in h are then solved 
for in a systematic perturbative fashion. We can also analyze the de Sitter case numerically; 
we find that we agree with the results of Basu and Vilenkin concerning the fact that for 
sufficiently large values of h, no wall solutions exist - when the horizon size is smaller than 
the domain wall width, a kink configuration cannot fit in the horizon. 

We next perform an expansion in the physical spatial coordinate that is nonperturbative 
in h, and use it to examine the behavior of the energy momentum tensor for the wall near 
the origin. This is of interest due to the possibility that a more natural form of infiation 
may occur in the core of the defect, where the field is trapped in the unstable vacuum [^. 

In curved spacetime, renormalization arguments suggest that, in general, a coupling to 
the Ricci scalar has to be introduced in the bare Lagrangian. We find that such coupling is 
responsible for very remarkable effects that result in a change of sign of the stress tensor at 
the origin. In fact, we see that this effect can also happen in fiat spacetime, if the so-called 
improved stress tensor is used. 

We derive the effective Lagrangian for the fiuctuations perpendicular to the wall. We 
find corrections beyond the "Nambu-Goto" action (the 3- volume swept by the worldsheet of 
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the wall) in a consistent expansion in derivatives. We also find non-trivial contributions to 
the surface tension from gravitational effects. 

We use this effective Lagrangian to address the question of whether the shape of the wall 
fluctuates, or whether it remains fiat at long distances and long times after the formation of 
the defect. Garriga and Vilenkin |^ have discussed this question to a certain extent, however 
we find that by using collective coordinates the long-wavelength fluctuations are identified 
with Goldstone modes and the instability found by them is avoided. 

In the next section we set up our ansatz and perturbation expansion. Section 3 deals 
with the energy momentum tensor for the wall in an FRW universe. We find that interesting 
behavior can occur by making use of the coupling ^TZcfP' of the field to the Ricci scalar TZ. 
In section 4, we analyze fluctuations of the wall. We differ from Garriga and Vilenkin in 
our treatment of translational invariance. The fact that the wall is located at a particular 
point along the 2;-axis, say, would appear to break this symmetry. However, translational 
invariance can be restored by use of the collective coordinate method. Doing this, we arrive 
at a true Nambu- Goldstone mode for the fluctuations, as opposed to one with negative mass^ 
as found by Garriga and Vilenkin. Several equal (comoving) time correlation functions that 
determine the behavior of the shape of the wall are computed, and we find that their spatial 
fall-off is much slower than in flat Minkowski spacetime at long times. This leads us to the 
conclusion that the surface remains "rough" at large separations. 

Section 5 contains our concluding remarks and some cosmological implications of our 
results. 

2 The Domain Wall 

Domain walls (interfaces) are field configurations with non-trivial topological boundary con- 
ditions. These boundary conditions demand that the field vary substantially within a typ- 
ical spatial scale that is usually determined by the (finite temperature) correlation length 
^ m^^, with m the, in general temperature-dependent, effective mass of the field. 

In FRW cosmologies there is another important scale, the microphysical horizon size 
dh = H^^ where H is the Hubble parameter. Causality implies that the region in which the 
scalar field can vary appreciably must be subhorizon sized, since only then can the boundary 
conditions that define the wall fit inside the horizon. Thus we expect that the notion of a 
domain wall will only make sense if H/m < 1. 

When H/m <^ 1, we expect that the domain wall profile will exhibit only small deviations 
from the Minkowski spacetime profile. We can then study the differences in a power series 
expansion in H/m = ^/dh- 
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Consider a scalar field with potential 

u = \{<p'-v^f + \n<j? (1) 

in a FRW spacetime with metric g^^ = diag(l, — a^, — a^, — a^), and Hubble parameter H = 
a/ a (where a = da/dT). The second term couples the field to the Ricci scalar TZ = 6{2H^ + 
H). Such a term is required for consistency as it will be generated by quantum corrections 
in general; the only fixed point is ^ = 1/6 which corresponds to conformal coupling [|l^. The 
equation of motion is 

where V is the derivative with respect to co-moving spatial coordinates {X,Y,Z}. It is 
convenient to define the following dimensionless quantities 

m = "\/Af , {x,y, z,t} = m{X,Y, Z,T}, h = H/m, 

r] = (f)/v, uj = a{t)z, f3 = 4 + H/H'^, fi"^ = 1 - Q^h\(3 - 2) . (3) 

Note that uj is a (dimensionless) physical coordinate. For power-law expansion (PLE) 
a{t) = t", we have j3 = A — l/n. In the radiation dominated (RD) and matter dominated 
(MD) fiat universes respectively, Pro = 2 and Pmd = 2.5. In de Sitter spacetime Pds = 4 
(or equivalently, we can take n = oo). 

We will assume a domain wall (kink) along the 2;-axis and that the corresponding field 
configuration is independent of the transverse coordinates, so that we are considering a fiat 
interface. 

The scale factor dividing the (comoving) Laplacian in eq. (H) suggests the following ansatz 
for the kink configuration: 

rj{z,t) = 'n{u;,h{t)). (4) 

Such a solution obeys 

(1 - h'u;'W' - (3h'u;i + /x^r/ - = \h'^ - -h'u;^ + (^--]h'^ (5) 



n"^ dh"^ n dh Km? nj dh 

where a prime means d/duj. The right-hand side (RHS) vanishes for de Sitter spacetime, 
and is C»(/i^) for PLE. 

In the case of de Sitter spacetime, the effect of non-zero conformal coupling ^ {i.e. ^ ^ 1) 
can be absorbed in a simple rescaling of variables: 

rj = fif], uj = (jj/n, h = jjih. (6) 

For PLE, this introduces additional Oiji^) terms on the RHS of eq. (^, since /i 7^ 0. 
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Except to discuss the energy- momentum tensor T^,y, we will hereafter work in the rescaled 
theory (and drop the bars), writing 

n > 2 2-, f, I ^ 3/0 (de Sitter) 

and taking "soliton boundary conditions" ■/^(ioo, /i) = ±1. All results for PLE are thus valid 
only to 0(/i2). 



2.1 Walls in de Sitter Spacetime: Numerical Solution 

The differential equation is, in general, a partial differential equation in two variables, 
and not easily amenable to numerical study. As we see in eq. (0), ignoring the RHS would 
introduce errors at 0{h^) for PLE. However, in de Sitter spacetime, where h is constant, 
we have an ordinary differential equation, for which we can find exact numerical solutions. 
Indeed, this was done by Basu and Vilenkin [Q, and our results appear to agree with theirs. 

Our results are plotted as thin curves in Fig. 1. The slope at the origin was chosen (by 
a shooting algorithm) by requiring smoothness through the point u = 1/h (marked with a 
tick). As found by Basu and Vilenkin 0, solutions only exist when < 1/4. (For PLE, if 
one insists on ignoring the RHS of eq. (^, solut ions only exist when h"^ < 1/(3.) 

We see the kink has been flattened out by the effects of the cosmological expansion. The 
slope at the origin ai = ri'{0, h) decreases with h"^ as shown in Fig. 2, vanishing at /i^ = 1/4. 
In other words, for a topologically stable solution to exist, the horizon size must be greater 
than the correlation length (by a factor of 2). 

2.2 Expansion in 

Although a numerical solution is available in the case of de Sitter, the general case is ex- 
tremely difficult to analyze numerically, since we now have a non-linear partial differential 
equation in two variables. However, when h"^ <^ 1//5, we can solve eq. (^ in a systematic 
perturbative expansion in h"^, and the first-order result will be valid for both de Sitter and 
PLE cosmologies. We write 

T]{co,h) = ri,ico) + J2h''^S^''\co) . (8) 

n 

Here 77^ is the kink configuration which solves eq. (|^) for /i = (flat spacetime): 

r]s{uj) = tanh(tu/V2) (9) 
This h = kink is shown as a dashed curve in Fig. 1. Note that //^(O) = 1/^2. 
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Since rjs has the correct asymptotic behavior, we must set 5*^")(±oo) = 0. To leading 
order in /i^, S{u;) obeys 



OsS^'^ = ^ + ^^'^ - = i(^) ^ vMu;' + m^)^ ■ (10) 



d5^ 

The solution may be found by elementary methods once a solution of the homogeneous 
equation is found. Such a solution is easily available, since translational invariance of the 
unperturbed differential equation guarantees that 

<5iH=r?:H = i=sech2(^^j, (11) 

is a "zero mode" of the operator Og of quadratic fluctuations around the kink conflguration. 

Another linearly independent solution 62{uj) having unit Wronskian with 5i{u!) can be 
found by elementary methods: 



52{UJ) = ^ 



sinh (V2uj) + 3tanh f ] + -^scch^ i —i= 



(12) 



where we have chosen this solution to vanish at the origin. Finally the flrst order solution 
can be written 



5( 



'\uj)=S2{uj) rS,{C)j{C)dC-S,{uj) rS2{C)j{C)dC + aS,{uj) + bS2{uj), (13) 
Jo Jo 

where the coefficients a, b have to be determined from the boundary conditions. Notice that 
Si{u!) is symmetric whereas j(a;), 62(00) are antisymmetric functions of cu. 

Note that ^2(0;) diverges exponentially as a; — > ±00. To render the solution finite and to 
satisfy the boundary conditions, we must choose: 

b=- Ts^iOjiOdC- (14) 
Jo 

The boundary condition at infinity does not determine the constant a. Notice that the 
term aSi{u;) represents a local translation of the wall. Clearly the freedom of choice of a 
refiects the underlying translational invariance. If we demand that the wall position be at 
the origin, we must demand that 5^-^^ (a; = 0) = 0; this leads to a = 0. 

Thus the final solution to this order that satisfies the boundary conditions and keeps the 
wall centered at the origin is given by 

S^'\uj) = -d2{uj) ^^{C)J{OdC-^,{uJ) 52(C)j(C)dC- (15) 

Ju) Jo 

The results for de Sitter spacctime (P = 4) are shown as thick curves in Fig. 1. S^^\uj) 
for /3 = 4 is plotted in Fig. 3. This approximation underestimates the distortion away from 
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the h = kink, and does not see the singularity at /i^ = 1/4, which is well beyond the regime 
of validity of the perturbative expansion in h^. 

It is straightforward to see that the structure of the perturbative expansion persists to 
all orders. In fact we find that 

with the same differential operator Os as in (p!OD, and where the source term obtained 
from the solutions up to order n — 1. Clearly the same functions 5i{uj) , 52(^^) generate the 
solutions, as these are the two linearly independent solutions of the homogeneous differential 
equation Os6 = 0. Thus the general form of the solution is the same as ([13|) with j^^^ 
replacing j and coefficients a^^\ fo*^"^ replacing a, b. 

By the same arguments presented before, a finite solution at infinity requires that 

roc 

= - / 6^{C)j^-\0dC. (17) 
Jo 

The source term constructed from iterations up to order (n — 1), and for n = 1 

it is given by derivatives of the unperturbed kink solution that vanish at u = ±oo. It is 
straightforward to prove by induction that the source term j(")(cj) is antisymmetric and 
vanishes exponentially as cu ^ ±oo. Therefore b^"''' is finite to all orders. Translational 
invariance suggests that a'-"^ = as before, and the the solution to order n is antisymmetric 
and this property ensures that the source term for the next order iteration is antisymmetric 
again. Thus we are led to the conclusion that the general form of the perturbative solution 
is given by eq. (^) with j{u) replaced by and the only complication in carrying out 

this program to any arbitrary order is to find the source term iteratively. Clearly this will 
be a highly non-local function because of the nested integrals, and the numerical evaluation 
will become more cumbersome in higher orders, but in principle this scheme will yield a 
consistent perturbative expansion. Clearly we have no way of determining the radius of 
convergence of such an expansion. 

2.3 Expansion in u 

The field configuration corresponding to a domain wall has most of the gradient and potential 
energy difference (with respect to the broken symmetry vacuum) localized in a spatial region 
of the order of the correlation length. Most of the contribution to the energy momentum 
tensor, after subtracting off the vacuum value, will arise from this small region around the 
position of the domain wall. This motivates us to obtain the kink profile near the origin in 
a power-series expansion in the coordinate u. This expansion is nonperturbative in h"^. 

The motivation for this expansion is the following: one would eventually like to study the 
full system of Einstein's equations and matter field in a self-consistent semiclassical manner. 



6 



In particular we have in mind the question of how the presence of domain walls affect the 
gravitational fields. Since for this task one only needs the behavior of the energy momentum 
tensor, and as we argued above most of the contribution to its components from a domain 
wall field configuration arise from a very localized region near the position of the interface, 
a power series expansion near the origin may capture most of the relevant physical features 
for an understanding of the back reaction of domain walls on the gravitational field. This 
argument becomes more relevant in view of the proposal of rapid inflation at the core of 
topological defects 

In the general case of de Sitter or power-law expansion we find the following behavior of 
the interface profile near the origin: 



6 V 120 20 / ^ 



The linear coefficient ai cannot be determined perturbatively, and will have to be found by 
solving the differential equation, its value being determined by requiring smoothness through 
uj = 1/h. We notice, however,that as /i^ — > from below, the first term dominates (even 
for large u) and the profile flattens out to a straight line. 

Since this expansion works best for /i^ close to 1/(3, but is correct only to 0{h?) for PLE, 
it is only useful in this case very close to the origin. 

Given a value for ai (from the numerical solution) as input, we display this series [through 
0{uj^y\ for the de Sitter spacetime {(3 = 4) case as the thin dot-dashed curve in Fig. 1. For 
h'^ ~ 0.2, where the solution is very flat, the small-a; expansion is valid out to well past 
ui = 1/h. We will use these results in the following section. 



3 The Energy- Momentum Tensor 

We next turn to an examination of the energy momentum tensor for the domain wall profile 
in an FRW spacetime. 

In general, the energy-momentum tensor of a scalar field configuration is given by 

r.. = ^^^^P^, C = \g-^{d^<P){d,<P)-U{ct>) (19) 



where g = det {g^j,^}, and we use 

sy^ _ 1 

Sg^"' 2 



-9 9,u (20) 



At this stage we restore the non-minimal coupling to the Ricci scalar, since, as we ar- 
gued above, such a term will be induced by renormalization if not present in the original 
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Lagrangian |]T0[. Its presence has two major effects: the first one is to modify the form of 
the energy momentum tensor, while the second is to modify the field configuration for the 
domain wall. Unlike the analysis of the equation of motion in the previous sections, for the 
case of the energy momentum tensor we can no longer rescale variables to get rid of ^, since 
{6n/5g'"') ^ 0. Let Uo be the potential for ^ = 0, so f/ = f/o + (^/2)7^02. Performing the 
variation of eq. we find: 



In our FRW metric, the non-zero components of = / {m?v'^) are: 



(21) 



i\2 



1 - (h^Y (,,') 



+ i[{2(3- b)h^rf + (2/? - 2)h'^ujr]r]' - 2 [l - (hu)^] [{r]'f + 7]r]"] } 



-1 
2 



+ ^ {(2/5 - 5)h^ri^ + {2(3 - 2)h'^ujT]i + 2{hujf [{r]'f + r^r/"] } 
-a\t)t^ = ait)huj{ivr-2^[ivr + VV"]} 



(22) 



As usual, for PLE there are corrections at 0{h^). 

The presence of the non-linear coupling to the Ricci scalar modifies the definition of 
the energy-momentum tensor. Even in the case of Minkowski spacetime, with h = 0, this 
modification of the energy-momentum tensor produces some peculiar behavior in near 
the origin, despite the fact that this coupling ^ does not affect the equations of motion. In 
this case, the wall solution is just rjs from eq. (|^), and the non-zero components of are: 



- - 2^ + ^cosh(V2cj) sech^(w/\/2) 



(23) 



For large enough ^, the energy density is negative at the origin: {t^){0) < for ^ > 1/2. For 
^ < 0, a region of negative energy density occurs away from the origin. 

We can substitute the small-tu expansion rj = aiU + 0{uj^) (which is correct even for 
7^ 0) into eq. (|22| ) to find the behavior at the origin. The non-zero components of t^{0) 
are: 



C(0) = t/(0) = t|(0) = ^ + |-2a?e, 
and the energy density at the origin is negative for 

, 1 1 
^>4 + 8^- 



(24) 



(25) 
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These results are consistent with eq. (p3D, where of = | [see equation (|])]. 

In Fig. 4 we show exact numerical calculations of for de Sitter spacetime with h"^ ^ 0. 
The behavior of for small values of h"^ (such as h"^ = 0.02 in Fig. 4a) is qualitatively similar 
to eq. (pSl). For larger h"^ (such as h"^ = 0.10 in Fig. 4b) another effect enters, namely that 
the wall profile fi{uj) is really determined hj h = h/ fi, as defined in eq. @. Wall solutions 
in de Sitter spacetime only exist for h"^ < 1/4, which implies 

< (26) 
4+12^ ^ ^ 

For example, with h"^ = 0.10, the wall profile flattens away to?7 = 0as^— >l/2 (where 
P ^ 1/4), and to°(^) = 1/4. 

The lesson that we learn from this analysis is that the coupling to the Ricci scalar can 
dramatically modify the behavior of the energy momentum tensor near the origin, and may 
be an interesting possibility for topological inflation at the core of defects [0]. These effects 
can even arise using the "improved" energy momentum tensor in Minkowski spacetime, 
if one starts with a general curved spacetime and then takes the flat limit. 



4 Fluctuations of the Wall 

Up to this point our study has focused on the description of a "flat" interface or domain 
wall, that is the fleld proflle varies only along the direction perpendicular to the domain wall 
but is constant on the perpendicular directions. However there will be fluctuations both 
quantum and thermal that will tend to distort locally the interface. An important question 
to address is the following: are these fluctuations "small" in the sense that the wall remains 
flat at long distances, or are the fluctuations important so that the wall becomes "rough"? 
We will assume that the system is at zero temperature and that only quantum fluctuations 
are important. 

We begin by deriving the effective action for the fluctuations of the interface and then 
proceed to calculate relevant correlation functions. 

We depart from the treatment of Garriga and Vilenkin [Q, in that we sacrifice explicit 
covariance to treat translational invariance in terms of collective coordinates; however, the 
final result for the action will be fully covariant. This procedure, borrowed from the usual 
scheme to quantize the collective coordinates associated with translations of soliton solutions 



ll| has many advantages. 

The position of the interface explicitly breaks translational invariance. However, a rigid 
translation of the interface should cost no energy due to the underlying translational symme- 
try. Thus the fluctuations perpendicular to the interface should be represented by massless 
degrees of freedom since, locally, they represent translations of the interface. These are the 
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capillary waves or Goldstone bosons |]T2|, |T^, |n| associated with the breakdown of the trans- 
lation symmetry. Translational symmetry is then restored by quantizing these fluctuations 
as collective coordinates [|n|. 

In terms of dimensionless comoving coordinates and the dimensionless field rj, the action 
is (d^a; is the comoving volume element in dimensionless units) 



/ = i/d^xa3(t) ifl^^ 1- 

XJ ^'\2\dtJ 2a?{t) 



1 



(27) 



We want to incorporate fluctuations of the interface solution (kink) and to obtain the 
effective action for the long-wavelength modes. In order to understand how to achieve this 
goal, it proves convenient to recall how this procedure works in Minkowski spacetime. 



4.1 Goldstone Modes in Flat Spacetime 



Consider first a static kink along the z-direction. As is well known |Tl| there is a zero mode 
of the linear fluctuation operator 

« (28) 

where Zq is the position of the kink and rjg is the kink solution. This mode corresponds 
to a translation of the position of the kink, because r]s{z — zq) + a Sot] ~ r]s{z — zq + a). 
Translational invariance guarantees that such a perturbation does not change the energy 
and thus is a zero mode of the linear fluctuation operator. Now consider a kink in three 
space dimensions, corresponding to a wall along the 2;- axis and centered at zq. A perturbation 
of the form 

6v,{x, y, z-z,) = ^^ ^^^(^-^0) ^.,-,..^, ^ ^29) 
with ■ x_L = Px^ + PyUi is an eigenmode of the linear fluctuation operator with eigenvalue 



vl + Vy [0 • These fluctuations of the interface correspond to the Goldstone bosons 



of the broken translational symmetry, and are the capillary waves of the interface [|T2i |T3i [14 |. 
The perturbed solution 

r]s{z - zq)+Y^ 5r]p{x, y,z- Zq) ^ ri,{z - f{x, y)) (30) 
p 

f{x,y) = Zo-Y.a,e^P^--'- 
p 

corresponds to a local translation of the interface. The correspond to "flat directions" 
in function space. Just as in the one- dimensional case, these Goldstone modes cannot be 
treated in perturbation theory, because arbitrarily large ap for pj_ — > can be accessed at 
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no cost in energy. We will borrow results from the one-dimensional procedure and quantize 



these modes as "collective coordinates" 15 



Besides these Goldstone modes there are massive modes corresponding to the higher 



energy states of the one- dimensional kink with dispersion relation E(pj_) = \/p± + Sm^ 



111 , [121 . Because of this gap in the energy spectrum, we can safely concentrate on the long- 
wavelength fluctuations of the interface and obtain an effective action for these Goldstone 
modes, treated as collective coordinates. The coordinates Op or the field f{x,y) are now 
fully quantized. In a path integral quantization procedure, collective coordinate quantization 
amounts to a functional integral over all configurations of f{x,y). Clearly this procedure 
restores translational invariance because now the field is invariant under z z+6, f{x, y) 
f{x, y) + S and the field f{x, y) is now functionally integrated with a translational invariant 
measure. Passing from the original integration variables in the functional integral to the new 
variables including f{x, y) involves a Jacobian which is seen to be unity to the order that we 
are working in (see below). For a thorough exposition of collective coordinate quantization 
in the path integral and Hamiltonian forms, the reader is referred to the original literature 



ri, ra, ra F 



4.2 Goldstone Modes in a FRW Cosmology 

After this digression in Minkowski spacetime we are ready to extend these observations to a 
FRW cosmology. Let ris{uj) be a solution to eq. (|^ (the /i-dependence will now be implicit), 
and let us take the small-gradient limit, in which 

fxx, fyy, ftt<^l; hf^, hfy, hft < 1 (31) 

{fxx ■ ■ ■ are second derivatives with respect to the dimensionless comoving coordinate x, 
etc.) This approximation is consistent with our purpose of studying the long- wavelength 
fluctuations on distances such that 1 <^ x, y ^ h^^. We look for a profile of the form 

V{z, f{x, y, t),t) = ris{^{z, X, y, t)) (32) 

Note that the substitution C,{z,x,y,t) = a(t)[z — f{x,y,t)] does not in general solve the 
equations of motion, but 

cf +\ f+\ [z-f{x,y,t)] 

t[z, X, y, t) = aft) , = (33) 

'I + fx + f' - a-{t)f! 



does, in our small-gradient approximation. 

The denominator in eq. (|33| ) has an important physical meaning. The function f{x,y,t) 
determines the position of the domain wall (interface). This function induces a metric g'^fj 
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on the 2 + 1 dimensional world-volume swept out by the wall, and we find, in terms of 
dimensionless comoving variables. 



VV^) = a2(t)y'l + /2 + /2_a2(t)/2 



(34) 



The effective action for the displacement field f{x,y,t) in the long- wavelength approxi- 
mation, is found by following the usual procedure |]T2|, |T^, which consists of computing 
the action for the profile (|32D with (^) (this is also identified with the effective action for 
the collective coordinate [|l^). After integrating by parts and discarding surface terms, we 
find 



+ 



Co 



d^x 



h h^' 

Co-C,{5 + 2-)- 



{l + f^ + f^-a'it)fi 



3/2 



(35) 



dt(l + f. + f'-a\t)fi 



a{ty 



a{ty 



where d^x represents the 2 + 1 dimensional comoving volume element, and we have used 
?7s(— = ^VsiO eliminate the integrations linear in C, as well as having defined: 



Co 



Co 




If 



This is the final form of the action for the fiuctuations of the interface. It exhibits the 
translational symmetry explicitly, and only contains derivative terms as is required of a 
Goldstone field. Co is identified with the (usual) fiat spacetime surface tension, and we see 
that curved spacetime effects induce a renormalization of this surface tension. The first term 
(proportional to \fg^) is recognized as the equivalent to the "Nambu-Goto" action, which is 
essentially the total "world-volume" associated with the fiuctuation field. The second term 
is thus a correction to the "Nambu-Goto" action; further corrections can be obtained in a 
systematic expansion in derivatives for the solution of the equations of motion. 

Expanding / to quadratic order in the fiuctuation field shows that the fiuctuations are 
massless: 



1 



/quad = Io + jld'x a\t) Co - C2(5 + 2h/h')h'/2 



1 



2a{t) 



-Xfl + fl) 



(36) 
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in contrast with the results found in reference 0. 

An important quantity that gives information about the behavior of the fluctuations of 
the interface is the vector normal to the interface: 



(37) 



Its correlation functions, to be computed below, will give information on whether departures 
from a fiat interface are significant. 

At this point we would like to compare our result with those of Garriga and Vilenkin. 
These authors found that the fluctuations perpendicular to the interface are associated with 
instabilities that manifest themselves as a tachyonic mass for these fluctuations. From the 
discussion above, we are led to conjecture that the appearance of this tachyonic mass is the 
result of a quantization that does not preserve translational invariance. 

Because our procedure does preserve translational invariance but is non-covariant in the 
intermediate steps, a direct comparison of our results is somewhat subtle. However, we 
can gain some insight by trying different parametrizations of our kink-profile. Consider the 
following the parametrization 

vizJix,y,t),t) = r]s{x{z,x,y,t)) (38) 



instead of that of equations (|3^j33| ). 

We see that x is not invariant under the the rigid translation F ^ F + Fq, z ^ z + Fq] 
instead a translation of the interface is compensated by a time dependent transformation F 
F + FQ/a{t), z ^ F + Jq. Because of the time derivative terms in the original action such a 
transformation is not an invariance of the action, which is then changed by terms proportional 
to h (time derivatives of the scale factor), thus breaking translational invariance. Following 
the same steps leading to the effective action found above, neglecting higher derivative terms, 
integrating by parts and rearranging terms we arrive at: 

Co r ^0 (3 + h/h' 



J d'x a\t){Co - ^2/^72)^1 + a-2(F2 + F^) - 



+ ^ I d'x F'' _ + 0{F''F,F,,). (40) 



Keeping only the quadratic terms in the action we find that the fluctuation field acquires 
a tachyonic mass, mjp = —Sh"^. This is the same value of the mass obtained by Garriga 
and Vilenkin Thus we conjecture that the scalar field that measures departures from a 
flat interface introduced by Garriga and Vilenkin is equivalent (at least to lowest order in 
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derivatives) to the scalar field F parametrizing the fiuctuations of the interface as in equa- 
tions (|38|j39|) . As explained above, this parametrization explicitly breaks rigid translational 
invariance in any FRW cosmology. The appearance of the mass term is understood as a 
consequence of this explicit breakdown of translational invariance, although this fact does 
not explain the tachyonic nature of the mass. 

Thus it seems to us that there are advantages and disadvantages in both formulations. 
Whereas the formulation of Garriga and Vilenkin is desirable in that it maintains explicit 
covariance, there is the feature of instabilities associated with the tachyonic mass of the fluc- 
tuations, which if our analysis is correct, indicates the breakdown of translational invariance 
in the quantization procedure. On the other hand, our formulation, in terms of collective 
coordinate quantization, sacrifices explicit covariance, although the final result is covariant, 
but explicitly treats translational invariance and its restoration via the collective coordinate 
quantization. The collective coordinates represent massless fields as a consequence of this 
translational invariance. 



4.3 Quantization of the Fluctuations 

Quantizing the fiuctuations of the interface allows us to answer some relevant questions about 
the dynamics of the interface. In particular we can answer the question that we posed at the 
beginning of the section that is whether the interface (wall) is fiat or strongly fiuctuating at 
long distances. In order to answer this question we must compute the correlation function 
of the vectors normal to the interface at long distances 1 ^ r < h^^ with r the distance on 
the two-dimensional surface of the interface. 

Since we are interested in long distance physics, we will only consider slowly varying 
fiuctuations of the fiuctuation field / and neglect higher derivative terms in the action, 
keeping only the quadratic terms (one can be brave and pursue a perturbative expansion 
but we will content ourselves here with a lowest order calculation) in the action. Repeating 
eq. m, 



1 

2A 



(41) 



The equations of motion for the field / are given by: 



fu + /t 4/1 - 



C2{^hh+h) 



Co - C2(5 + 2h/h^)hy2 



A/ 



(42) 



with A being the two-dimensional Laplacian. 

In the most general case, the time dependence of the above equation is far too complicated 
to pursue analytically and one would have to resort to numerical integrations. Thus we 
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concentrate on the case of de Sitter expansion with scale factor a{t) = e^^ that allows an 
analytic treatment. 

The fluctuation field is expanded in terms of creation and annihilation operators and the 
mode functions which are solutions of the above equation of motion: 



f{x,t) 



a 



1 
A 



Cn ^ SCo 



(43) 
(44) 



where the "surface tension" a has been absorbed in the definition of / to make it canonical 
and A is the (comoving) area of the (planar) wall. 

Performing the change of variables on the mode functions Vp = e~^'^*Xp(t) the equation 
for Xp reads: 



- ^h' = 0. 



(45) 



whose solutions are linear combinations of the Bessel functions. The mode functions Vp are 



Vp{t) 




-2ht 



Ap H2 



(1) / 



h 



Bp H2 



(2) / 



h 



(46) 



where the coefficients Ap, Bp are arbitrary so far. 

Imposing canonical commutation relations between / and its canonical conjugate mo- 
mentum Uf = 6C/6ft = crftd^ leads to the relation 



\A^\ 
\^p\ 



\Bl 



1. 



(47) 



It is a well known feature of quantization in curved spacetimes that a choice of Ap 
corresponds to a choice of vacuum state. Although we do not have a physical criterion to 
pick a particular vacuum state, we will choose the Bunch-Davies vacuum for simplicity. 
Such a choice implies 

(48) 



Bp = 0. 



Without loss of generality we can take Ap = 1. Finally, the field / is expanded in creation 
and annihilation operators with respect to the Bunch-Davies vacuum state as: 



-2ht 



TT 



AhaA 



E 



-ht' 



h 



-ht' 



h 



(49) 



As mentioned above, we are interested on the long-wavelength fiuctuations of the interface 
at long distances and at long times after its formation. Thus we will study the regime 
t ^ h^^ ; hr < 1. Furthermore there are physical cutoffs that we must introduce: the 
(comoving) wavelengths cannot be bigger than the horizon, and because the nature of our 
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approximation cannot be shorter than the correlation length (the expansion is in derivatives, 
thus valid for slowly varying fields on the scale of the correlation length) . Therefore integrals 
over wavevectors will be restricted to the interval h < p < 1. 

The fiuctuation field f{x,y,t) measures the departure from a fiat interface. Thus a 
quantity of interest is the correlation function of the vector normal to the space-like interface. 
This normal vector is obtained from the induced metric on the two-dimensional surface and 
given by 

n = ^"^^7/" ^2 - (-/., -fy, 1 - fj2 - fy/2). (50) 

Vl + /' + fy 

The equal-time two-point correlation function is given by 



< n{x,t) ■ n{y,t) >~ 1 + dpp^[Jo{pr) - 1] 



(1) ^P^_^ 



(51) 



notice from eq. (|50| ) that < n{x, t) ■ n{x, t) >= 1 (to the order considered) and that the 
correlation function does not require any short distance subtraction. Its long time, long 
distance behavior is found to be 

<n{x,t) ■n{y,t) lnr/2H (52) 

a 

where ■ ■ ■ stand for sub-leading terms that fall off fast at large r. 

It is instructive to compare this result with that in Minkowski spacetime, for which 
< n-n >^ l — r~^. Clearly the interface is "rougher" in de Sitter spacetime. The result ( ^2]) 
raises the very interesting possibility of anomalous exponents in the correlation function. 
Notice that for large r the logarithmic (infrared) singularities become strong and eventually 
would have to be resummed in order to obtain meaningful long distance behavior. Assuming 
such a resummation of the lowest order result we obtain the long distance behavior: 

< n{x,t) ■ n{y,t) > ^ r"" (53) 

a = . (54) 

a 

This result has interesting implications. In particular, we find that at long times and 
distances the vectors normal to the interface are uncorrelated and the interface is "rough" 
rather than approximately fiat with small fiuctuations that fall-off rapidly at long distances. 

This situation is very similar to that of the X-Y model and other models in statistical 
mechanics in 1 + 1 Euclidean dimensions (typically massless scalar field theories) in which 



logarithmic singularities sum up to anomalous dimensions ||20| much in the same way. This 
would indeed be a tantahzing possibility that needs to be studied further, perhaps by keep- 
ing higher order terms in the effective action and resumming using renormalization group 
arguments. This is certainly beyond the scope of this article. 
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Another quantity of interest that measures properties of the interface is the distortion 

i: 

D\x, y- 1) =< [/(x , t) - f{y, t)f >= 2 < f{x, t) > -2 < fix, t)f{y, t) > (55) 

We evaluated this correlation function for large times and distances {t 1/h ; r^l^ hr) 
and found 

D\X,Y) ^ li^r^lnr. (56) 
cr 

This result is consistent with that of the normal- normal correlation function (|5^) . 

This correlation function shows once again that the interface is strongly fluctuating at 
long distances and cannot be considered flat. 



5 Conclusions 

In this article we have focused on the properties of domain walls in general FRW Universes. 
These defects will appear during phase transitions in typical inflaton theories and ultimately 
drive the dynamics of the process of phase separation. They may also contribute to density 
fluctuations and perhaps to structure formation. 

Our results present a rich picture of the properties of these defects. After analyzing 
numerically the case of de Sitter spacetime and assuming that the boundary conditions on 
the field configuration require that the horizon size be larger than the typical correlation 
length we set up a systematic perturbative expansion in powers oi h = H/m and worked 
out in detail the first order correction from curved spacetime effects. We see then that 
even in situations in which the Universe is expanding slowly enough to allow us to set 
up our perturbation theory for defects, new phenomena occurs with definite cosmological 
implications. 

In particular, we argue that in a general FRW spacetime, consistency (renormalizability) 
of the theory requires a coupling to the Ricci scalar. With this coupling the energy momen- 
tum tensor of defects can show some remarkable new features, including a negative energy 
density at the origin. Clearly this observation brings interesting possibilities for topological 
inflation near the center of the defect that must be studied within the full set of Einstein's 
equations in the presence of this configuration. For some values of the coupling to the Ricci 
scalar, we find that the energy density at the center is lower than at a distance of the order 
of the correlation length away from the center. 

Such a behavior of the energy momentum tensor may provide an interesting mechanism 
for scalar density perturbations. 

Using collective coordinate quantization of the fluctuations perpendicular to the wall, we 
obtained the effective Lagrangian for these fluctuations in the long-wavelength approxima- 
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tion. There are several noteworthy features of this effective Lagrangian: curved spacetime 
effects renormahze the surface tension, and we find systematic corrections to the "Nambu" 
action. 

We have also shown how to quantize the fluctuations about a completely flat wall, and 
how the collective coordinate procedure yields the expected Goldstone mode associated with 
translational invariance. This analysis then allows us to compute correlation functions and to 
see that, at least in the de Sitter case, fluctuations are fairly strong at long distances and long 
times after formation of the defect leading to conclusion that the surface tends to be rougher 
than in flat spacetime. These results and the method developed for the quantization of the 
fluctuations may yield some insight on novel mechanisms to study density perturbations. 

What we have done for walls can also be done for other defects, such as cosmic strings. 
These should exhibit some interesting behavior, especially for the fluctuations, since there will 
now be collective coordinates associated with the internal U{1) symmetry whose spontaneous 
breaking gives rise to the string. Interesting possibilities are corrections to the string tension 
from curved spacetime effects and also corrections to the "Nambu-Goto" action from higher 
derivative terms, and perhaps novel behavior of the energy momentum tensor at the core of 
the defect when the coupling to the Ricci scalar is introduced. We are currently studying 
these and other issues. 
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Fig. 1: Domain wall profiles for de Sitter space, h? = {■!, -2, .24}. Dash=kink, thin=exact, 
thick=sniall-h^ , dot-dash.=small-u} . uj = ±l/h is marked by vertical ticks. 
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Fig. 3: for de Sitter space. 
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